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PROGRAM SUMMARY 

Title of program: Sphinx - HEP MONTE CARLO 
Catalogue number: 

Program available from: Oliver Martin, 

|http : / /th . physik . uni-f rankf urt . de/h ^ mart in/ mart in/ sphinx . html 



Computer: IBM RS6000; PentiumPro based PC, other machines with FORTRAN 77 
and sufficient capacity should be able to run this program 

Installation: Institut fur Theoretische Physik, Goethe-Universitat, Frankfurt 
am Main, Germany 

Operating system: tested under AIX and LINUX but does not depend on the 
particular operating system 

Programing language used: FORTRAN 77 

High speed storage required: 5 Mbytes 

Number of bits per byte: 8 

Number of lines in combined test deck: 19296 

CPC subprograms used: Jetset 7.3 

Keywords: polarized nucleon-nucleon scattering, high energy physics, Monte Carlo 
simulation 

Nature of physical problem: This program can be used to simulate polarized 
nucleon - nucleon collisions at high energies. Spins of colliding particles are 
taken into account. The program allows to calculate cross sections for various 
processes. 

Method of solution: The existing Monte Carlo program Pythia 5.6 has been 
modified to incorporate spin effects. The program incorporates all features 
of Pythia. 

Restrictions on the complexity of the problem: Spins of colliding nucleons must 
be parallel to the collision axis. 

Typical running time: ~ 0.01 s on PentiumPro 200 based PC, depends strongly 
on kinematical cuts 
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References: Reader has to refer to Pythia manual M for additional details. 
Further informations can be obtained from: Oliver Martin, Institut fur The- 
oretische Physik, Goethe Universitat, Robert-Mayer-Str. 8-10, 60054 Frank- 
furt am Main, Germany, e-mail: martin@th.physik.uni-frankfurt.de . 
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LONG WRITE-UP 



1 Introduction 

The detailed investigation of hadronic spin effects, most notably in deep- 
inelastic lepton-nucleon scattering, has proven to allow for extremly sensi- 
tive tests of QCD. One of the most interesting issues is the possible con- 
tribution of the axial-vector anomaly. This anomaly plays a central role in 
modern quantum-field-theory and spin effects might offer the only possibil- 
ity to actually detect it. To estimate its contribution one has to know the 
distribution function for polarized gluons Ag(x,Q 2 ). It is argued that the 
total spin carried by gluons could be large, because a s (Q 2 )Ag(Q 2 ), with 
^9{Q 2 ) — lo Ag(x, Q 2 ), is a renormalisation group invariant, such that if 
Ag(Q 2 ) is a few percent at the constituent quark level it could be very size- 
able at large Q 2 . In deep- inelastic scattering it is impossible to distinguish 
between the regular quark contribution and the gluon contribution via the 
anomaly. There is actually no conceptually sound way to separate them. In 
exclusive or semi-inclusive lepton-hadron scattering channels it is in princi- 
ple possible to discern a gluonic and a quark contribution, but the by far 
most direct and conceptually clean way are polarized proton-proton colli- 
sions. RHIC [0] offers the option to study polarized proton collisions up to 
y/s = 300 GeV which would be ideal for this purpose. To analyse this option 
in detail, to plan the experiments, and to determine the reachable precission 
a detailed Monte Carlo program is needed, which we provide with this pub- 
lication. This program can handle, however, only longitudinal polarization. 
All polarization effects in proton-proton reactions share the same problem. 
As only a small fraction of the partons is polarized the signal to background 
ratios are typically of the order percent, such that a detailed simulation is 
necessary to extract quantitative results. 

Because a state-of-art Monte Carlo program is usually the result of many 
years of development we modified an existing, widly used program called 
Pyhtia |]J . We have been interested only in those parts of the program that 
correspond to nucleon-nucleon collisions. Consequently only these routines 
were changed. The resulting polarized version of Pythia is named Sphinx, 
an acronym for Simulator of Polarised Hadronic YNteract(K)ions. 

It is interesting to note that Sphinx is able to function in unpolarized 
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mode as well. In that case it is functionally equivalent to Pythia. This 
means that results will be the same but because some data structures inside 
the program have been modified the program requires more space. 

The main problem was to introduce the changes in a way compatible 
with the basic structure of Pythia, e. g. the way in which Pythia evaluates 
the cross sections or simulates the effects of higher order QCD corrections. 
Pythia is based on the parton model. The main assumption is that the 
nucleon may be represented by an incoherent mixture of quarks and gluons 
carrying portions of nucleon momentum. Their momentum distributions 
are described by parton distributions. In the parton model the nucleon is 
composed of quarks and antiquarks of various flavours and of gluons. In the 
polarized case each parton comes in one of two helicities: righthanded and 
lefthanded. Therefore the number of parton components has to be doubled, 
corresponding to the two degrees of freedom for longitudinal polarization. 
The most natural idea was to generalize the notion of flavour to account 
for the polarization. So instead of d-quark we have lefthanded d-quark and 
righthanded d-quark etc. Because in Pythia gluons are treated as a special 
"flavour" the scheme applies to them as well. As far as the general structure 
of the program is concerned the only complication of the polarized code is 
that we have more "flavours" and this is rather easy to implement. 

Polarization is followed in Sphinx only up to the hard partonic inter- 
action, i. e. the hadronic cross section and the higher order corrections in 
the initial state shower are evaluated spindependently whereas all final state 
interactions as fragmentation and decays are treated spinindependently. No 
theoretical model for polarized fragmentation exists. It is, however, generally 
accepted that it should not depend on longitudinal polarization. For trans- 
verse polarization the situation would, however, be different. One knows 
experimentally that substantial effects exist, e. g. a strong coupling between 
transverse spin and transverse momentum. 

The polarization effects have been implemented in the following parts of 
the program: 

1. Parton distributions. At the moment we supply seven sets of po- 
larized parton distributions: Altarelli-Stirling [BJ, and six sets from 
Gehrmann&Stirling ||. First order Altarelli-Parisi [|| evolution is 
taken into account in all cases. 

2. Hard processes. The processes currently implemented in the polarized 
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mode are summarized in Table [I]. Other processes may be used as well 
but the results will be averaged over polarizations. The cross sections 
for all polarized processes are given in the literature j7J, and we checked 
them. 

3. Initial state showering. The polarized case was not discussed in this 
context before. We obtained all necessary formulae. 

4. Documentation. The event listing has been modified such that it dis- 
plays polarization information about the interacting particles up to the 
hard partonic interaction. 

Hereafter we do not try to explain the structure of Pythia 5.6 and 
we indicate only the most significant modifications which lead to Sphinx. 
Readers unfamiliar with Pythia should first consult the appropriate manual 

I- 

2 From Pythia to Sphinx: modifications 

Our program is an extension of Pythia version 5.6. The Pythia 5.6 pro- 
gram heavily relies on Jetset 7.3 subroutines H. These subroutines are 
mainly used to simulate processes that take place after the hard partonic 
interaction, i.e. fragmentation and decays, or provide spinindependent ma- 
nipulations as for example Lorentz-transformations. Because polarization 
effects in the final state are neglected in Sphinx, the Jetset program could 
be left untouched. Only the set-up of an event listing which is also a task of 
Jetset in a Pythia run has been taken over from Sphinx in the polarized 
mode. The event listing in SPHINX provides then in addition information 
about the polarization flow. For that purpose two new subroutines have 
been added which are modifications of the corresponding Jetset routines 
(see below). By this means Sphinx provides a suitable interface to Jet- 
set such that the programs may be linked without any modifications in the 
latter. There is only one exception: the Jetset subroutine LUGIVE, which 
returns the values of all Jetset and Pythia common block variables cannot 
be longer used, because in Sphinx the dimensions of some Pythia arrays 
have been enlarged so that the formats no longer match. 
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Table 1: List of processes implemented in the polarized mode. 



I SUB 


Process 


Comment 


1 




quark-antiquark annihilation into virtual 

7*/Z° 


2 


m - 


annihilation into charged vector boson 


11 


QaQa — > QaQa 


(anti-) quark - (anti-) quark scattering; anni- 
hilation diagram not included 


12 


q^i -> ftft 


annihilation process 


13 


ftft -»• 


annihilation into gluon pair 


14 


ftft -»• #7 


annihilation into gluon and prompt 7 


15 


gift -> 


annihilation into gluon and Z° 


16 


gift -> fi'W^ ± 


annihilation into gluon and W ± 


18 


gift -»• 77 


annihilation into 7 pair 


19 


gift -> 7Z 


annihilation into 7 and Z° 


20 


gift -> 7^ /± 


annihilation into 7 and W 


28 


Qi9 -> ft# 


(anti-) quark - gluon scattering 


29 


ft£ -> ft7 


prompt 7 production in (anti-) quark - gluon 
scattering 


30 


$3 -> ftZ 


Z° production in (anti-) quark gluon 
scattering 


31 


ft^ -> ft-w^ 


W ± production in (anti-) quark gluon 
scattering 


53 


-> gfcft 


gluon fusion 


68 


99 -> 


gluon - gluon scattering 
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In many places inside Pythia partonic data are stored in arrays indexed 
by flavour. Inclusion of helicities means doubling the number of "flavours" 
as described in Section |I[ For that reason we had to enlarge the arrays. The 
unpolarized mode has been preserved throughout the program and works 
exactly as the original code. All modifications are implemented in a treelike 
structure with several branching points where the user may choose between 
the polarized and unpolarized modes. The default is always unpolarized. An 
interesting point is that one can combine polarized and unpolarized treat- 
ments for different processes and for various stages of the same process. This 
has to be done with care and one has to be aware of this not always being 
justified from the physical point of view. The branching points were imple- 
mented by means of the IF-ELSE structure which directs the program flow 
according to a specific switch. Each functional part (usually a subroutine) 
has its own local variable called IPOL. This variable controls the mode and 
its value depends on some parameters^ MSTP(x) and NSUB(x) to be defined 
by the user in the main program. For example, it is possible to run the code 
with polarized partonic cross section and polarized parton distributions but 
unpolarized initial state showering etc. This solution provides flexibility but 
has to be used carefully. The major disadvantage of IF-ELSE constructions 
is that parts of the code are multiplied. This makes servicing of the code 
more cumbersome because one has to introduce the same change in different 
places at once. In addition the resulting code became quite lengthy. 

Sphinx is as Pythia a "slave-system", i. e. it consists only of callable 
subroutines where two of them (PYINIT and PYEVNT) have to be called by 
the user to perform the event generation and a few other, e. g. DPLIST, 
PYSTAT, etc., could be called to obtain further event information. The rest 
of the subroutines is of internal use. Therefore the user has to supply a main 
program where all relevant parameters and switches have to be specified and 
these subroutines have to be called. Because this structure is the same as 
in Pythia the reader is referred to |1J again for details about the general 
Pythia parameters. In this article only the new parameters in Sphinx 
are discussed. The modifications made in the individual subroutines are 
described in Section |[ In Table [I] we list all new parameters introduced to 
control the polarized mode. Finally in Section |] examples for a main program 
and the result of the corresponding test runs are presented. 

^ee Tabic g 
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Table 2: Parameters controlling the polarized mode 



Parameter 


Description 


Default 


MSTPC171) 


=0 
= 1 
=2 


beam polarization 
unpolarized 

polarization in +z direction 
polarization in — z direction 





MSTP(172) 


=0 
= 1 
=2 


target polarization 
unpolarized 

polarization in +z direction 
polarization in — z direction 





MSTPC175) 


use of polarized parton distribution in polar- 
ized initial state shower (MSTP(176)=1) 

=0: unpolarized distribution; for testing only, do 
not use! 

= 1: polarized distribution 


1 


MSTP(176) 


initial state showering mode 
=0: unpolarized 
= 1: polarized 
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Table 2: Parameters controlling the polarized mode, continued 



Parameter 


Description 


Default 


MSTP(177) 


set of polarized parton distributions Aq and 

Ag used; in addition one has to specify an 

unpolarized set as in standard Pythia 
=f) - An = D and An = D (no polarization ) 
= 1: fake polarization, built up from unpolarized 

distribution according to Ag = MST ^ 78 ^ q 
=2: Altarelli-Stirling parametrization 0; data 

file altsti.dat required 
=3: Gehrmann- Stirling parametrization 

lo set a; data file gestloa.dat required 
=4: Gehrmann-Stirling parametrization || 

lo set b; data file gestlob.dat required 
=5: Gehrmann-Stirling parametrization 

lo set c; data file gestloc.dat required 
=6: Gehrmann-Stirling parametrization nlo 

set a; data file gestnloa.dat required 
=7: Gehrmann-Stirling parametrization || nlo 

set b; data file gestnlob.dat required 
=8: Gehrmann-Stirling parametrization || nlo 

set c; data file gestnloc.dat required 
=9: fake polarization, for testing only, do not use! 
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Table 2: Parameters controlling the polarized mode, continued 



Parameter 



Description 



MSTPC178) 



percentage of fake polarization for 
MSTP(177)=1 



MSTPC180) 



mode selection (master switch) 
=0: unpolarized mode; this value overrides all 

other polarization switches 
=1: polarized mode 



NSUB(ISUB) 



mode for subprocess ISUB 
=0: unpolarized treatment 
=1: polarized treatment 



NSEL 



menue of polarized processes 



=1 


ISUB = 


1 1 , 12 , 13 , 28 , 53 , 68 switched on 


=10 


ISUB = 


14, 18,29 switched on 


=11 


ISUB = 


1 switched on 


=12 


ISUB = 


2 switched on 


= 13 


ISUB = 


15,30 switched on 


=14 


ISUB = 


16,31 switched on 
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Table 3: Internal variables storing polarization information 



Variable 


Description 


Com. Block 


MINTC311) 


=0 
=1 
=2 


beam helicity 
unpolarized 
positive helicity 
negative helicity 


PYINT1 


MINTC312) 


=0 
=1 
=2 


target helicity 
unpolarized 
positive helicity 
negative helicity 


PYINT1 


MINTC313) 


=0 
= 1 
=2 


helicity of shower initiator on beam side 
unpolarized 
positive helicity 
negative helicity 


PYINT1 


MINTC314) 


=0 
= 1 
=2 


helicity of shower initiator on target side 
unpolarized 
positive helicity 
negative helicity 


PYINT1 


MINTC315) 


=0 
= 1 
=2 


helicity of hard interacting parton on beam 
side 

unpolarized 
positive helicity 
negative helicity 


PYINT1 


MINTC316) 


=0 
=1 
=2 


helicity of hard interacting parton on target 
side 

unpolarized 
positive helicity 
negative helicity 


PYINT1 


MSTP(179) 


switch off polarization temporarely in 

PYSIGH and PYSTFU resp. 
=0: no action 
—1: switch off polarization 


PYPARS 
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Table 3: Internal variables storing polarization information, continued 



Variable 



Description 



Com. Block 



ISIGH(1000,6) 

ISIGH(I,1) 

ISIGH(I,2) 

ISIGH(I,3) 

ISIGH(I,4) 

ISIGH(I,5) 

ISIGH(I,6) 



hard scattering information of Ith line 
particle code of Ith line on beam side 
particle code of Ith line on target side 
colour flow 

helicity of Ith line on beam side 
helicity of Ith line on target side 
not used 



PYINT3 



KD(I) 



polarization/helicity of Ith line 
no polarization/helicity 
positive polarization/helicity 
negative polarization/helicity 



DPYPOL 



XSFX(2, -40:40, 0:2) 



XSFX(JT,KFL,0) 
XSFX(JT,KFL,1) 
XSFX(JT,KFL,2) 



x times parton distribution for given 
x and Q 2 of flavour KFL = -40:40 
and helicity KFLD = : 2 on beam side 
(JT=1) and target side (JT=2) resp. 
unpolarized 
positive helicity 
negative helicity 



PYINT3 
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Table 4: Parameters controlling the choice of the unpolarized parton distri- 
bution functions 



Parameter 


Description 


Default 


MSTP(51) 




Qpf of nirnol £\ T*T7Pri ni^Tfor) nKtrihntiriTici" cpp 


12 






cLlol) FlDl.rV.OZ J 






— i 

X 


F,fTT,01 cjpf MQRfi nnrtafprl Afpr^inn T,D1 






=2 


EHL02 set (1986 undated version LO) 






=3 


Duke-Owens set 1 (LO) 






=4 


Duke-Owens set 2 (LO) 






=5 


CTEQ2M (NLO) 






=6 


CTEQ2MS (NLO) 






=7 


CTEQ2MF (NLO) 






=8 


CTEQ2ML (NLO) 






=9 


CTEQ2L (LO) 






= 10 


CTEQ2D (NLO) 






= 11 


Gliick, Reya, Vogt (LO) 






= 12 


Gliick Reva Voet 94 (LO) 






= 13 


Martin, Roberts, Stirling MRS(AP) (NLO) 






= 14 


Martin, Roberts, Stirling MRS(Rl) (NLO), 








has not been tested yet! 






= 15 


Martin, Roberts, Stirling MRS(R2) (NLO), 








has not been tested yet! 




MSTP(52) 


= 1 
=2 


choice of proton parton-distribution-library 
the internal Sphinx one, with parton distri- 
butions according to MSTP (51) above 
the Pdflib one, with the Pdflib (version 4 
and higher) NGROUP and NSET numbers to 
be given as MSTP(51)=1000*NGR0UP+NSET. 
In order to make use of this option PDFLIB 
must be linked and the subroutines STRUCTM 
and PDFSET at the end of the program should 
be commented out. 


1 
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3 Common Blocks and Subroutines 



In the following the modified and the new subroutines that have been specif- 
ically created for Sphinx are described in more detail. The general tasks of 
the subroutines themselves as well as the unchanged parameters and variables 
are not explained, because they are the same as in Pythia. The reader is 
asked again if occasion arises to consult |I| to obtain the needed information. 
We restrict our explanation to the new aspects in Sphinx. The purpose of 
the modifications is indicated and the new parameters, switches, and internal 
variables are listed. Meaning and possible values of the new parameters are 
given in Table [I]. To incorporate polarization the following common blocks 
have been enlarged and replace the corresponding Pythia common blocks 
or are added: 

• C0MM0N/PYINT3/XSFX(2, -40: 40,0: 2) ,ISIG( 1000,6) ,SIGH(1000) 

• C0MM0N/PYSUBS/MSEL,NSEL,MSUB(200) ,NSUB(200) , KFIN (2 , -40 : 40) , 
CKIN(200) 

• C0MM0N/DPYP0L/KD(4000) 

Information about the new internal variables and enlarged arrays can be 
found in Table |l|. In addition it is shown how the local polarization switch 
IP0L is built up in the different subroutines. Only the polarized case (IP0L=1) 
will be discussed, because in the unpolarized case (IP0L=0) each subroutine 
works exactly as the corresponding Pythia subroutine. The not mentioned 
subroutines of Sphinx are the same as in Pythia. 

MAIN PROGRAM 

Purpose: to set up the polarized event generation. The variables which 
have to been set are listed in Table |l|. 

Remarks: Examples of a main program are given in Section |j. 
SUBROUTINE PYINIT (FRAME , BEAM , TARGET , WIN) 

Purpose: to display Sphinx header; to check partially the availability of 
the desired polarization scenario, i. e. to control that the master switch for 
polarization MSTP(180) is set properly, the selected partonic subprocesses 
can be treated polarized and to control and compose the polarization menue 
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via NSEL; to call DPLIST instead of LULIST (see below). 
New parameters: MSTP(180), NSEL, NSUB(ISUB) 
Internal Polarization switch: IP0L=MSTP(180) 

Remarks: The interface to the CERN parton distribution library was up- 
dated and works now with PDFLIB 4 and higher versions. Furthermore, if a 
not allowed scenario has been chosen, the programs stops with an appropri- 
ate error message. The PAKPDF parton distribution library is not supported 
anymore. 

SUBROUTINE PYEVNT 

Purpose: to start polarized event generation; to call DPEDIT instead of 

LUEDIT (see below). 

New parameters: MSTP(180) 

Internal polarization switch: IP0L=MSTP(180) 

SUBROUTINE PYINKI (CHFRAM , CHBEAM , CHTARG , WIN) 

Purpose: to check availability of the desired hadronic polarization scenario, 
i. e. to control that the selected hadron can be treated polarized and to 
verify that the polarization is longitudinal; to store the polarization of beam 
and target for the event listing in KD(1) and KD(2) and for internal use in 
MINT(311) and MINTC312). 

New parameters: KD(I), MSTP(171), MSTP(172), MSTP(180) 
New internal variables: MINT(311), MINT(312) 
Internal polarization switch: IP0L=MSTP(180) 

Remarks: At present only nucleons and hyperons and their antiparticles 
can be treated polarized. If a not allowed scenario has been chosen, the 
programs stops with an appropriate error message. 

SUBROUTINE PYRAND 

Purpose: to adapt PYRAND to the new environment - all relevant arrays 
which have been enlarged or added to the common blocks in other subrou- 
tines are modified here as well; to extend event shape selection to incorporate 
helicities; to store helicities of the partons entering the hard interaction ac- 
cording to 

• MINT (313): helicity of the beam parton for use in the initial state 
showering subroutine; 
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• MINT (314): helicity of the target parton for use in the initial state 
showering subroutine; 

• MINT (315): helicity of the beam parton; 

• MINT (316): helicity of the target parton. 

New parameters: MSTP(180), NSUB(ISUB) 

New internal variables: MINT(313), MINT(314), MINT(315), MINT(316) 
Internal polarization switch: IP0L=MSTP(180) xNSUB(ISUB) 
Remarks: Note that MINT(313)=MINT(315) and MINT(314)=MINT(316) 
but the values of MINT (313) and MINT (314) are changed later by the ini- 
tial state shower in PYSSPA. 

SUBROUTINE PYSCAT 

Purpose: to adopt PYSCAT to the new environment (see PYRAND); to store 

helicities of the partons entering the hard interaction; to fill lines 1, 2 and 5, 

6 in the event listing with polarization information (see below). 

New parameters: KD(I), MSTP(180), NSUB(ISUB) 

New internal variables: MINT(315), MINT(316) 

Internal polarization switch: IP0L=MSTP(180) xNSUB(ISUB) 

SUBROUTINE PYSSPA ( IPU1 , IPU2) 

Purpose: to perform polarized initial state showering, helicity dependent 
GLAP evolution equations are used in the backward evolution algorithmus; 
to enlarge all relevant array in an appropriate manner to incorporate polar- 
ization; to check proper selection of the polarized initial state shower scenario, 
i. e. to control that the MSTP(175) and MSTP(176) are set correctly; to store 
the helicities of initial state shower initiators (MINT (313), MINT (314)). 
New parameters: KD(I), MSTP(171), MSTP(172), MSTP(175, MSTP(176), 
MSTP(180), NSUB(ISUB) 

New internal variables: MINT(313), MINT(314), MINT(315), MINT(316) 
Internal polarization switch: IP0L=MSTP(180) xNSUB(ISUB) xMSTP(176) 
Remarks: At the present stage only QCD shower can be treated polarized, 
QED showering has to be done in the unpolarized manner. The combination 
MSTP(175)=0 and MSTP(176)=1 allows to simulate polarized showering with 
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the use of unpolarized parton distributions. This option is just for testing 
and should not be selected by the user! If MSTP(175) or MSTP(176) are set 
improperly, the programs stops with an appropriate error message. The in- 
ternal variables MINT (313) and MINT (314) are changed to their final values 
in this subroutine. 

SUBROUTINE PYMULT(MMUL) 

Purpose: to switch off polarization in PYSIGH (set MSTP(179)=1 temporally) 
when called from PYMULT even in a polarized run, because multiple interaction 
cannot be treated polarized at the moment. 
New parameters: MSTP(179) 

SUBROUTINE PYREMN(IPU1 , IPU2) 

Purpose: to adopt PYREMN to the new environment (see PYRAND); to fill lines 
3, 4 in the event listing with polarization information (see below). 
New parameters: KD(I) 

SUBROUTINE PYSIGH 

Purpose: to evaluate the helicity dependent hadronic cross sections by con- 
volution of the helicity dependent parton distributions with the helicity de- 
pendent partonic cross sections; to supply the subroutine with the helicity 
dependent partonic cross sections. 

New parameters: MSTP(171), MSTP(172), MSTP(179), MSTP(180), 
NSUB(ISUB) 

Internal Polarization switch: 

IP0L=MSTP ( 180) x NSUB (ISUB) x ( 1-MSTP ( 179) ) 

Remarks: PYSIGH will always run in the unpolarized mode when it is called 
by PYMULT which sets MSTP(179)=1 in PYSIGH temporally. Evaluating the 
spindependent hadronic cross sections one has to notice that the hadrons are 
specified according to the spin, whereas the partons are labelled by their helic- 
ities. The spin is defined relative to the collision axis and the beam is assumed 
to move in the positive direction. For that reason the helicities at the target 
side are opposite to the polarizations. Hence the target labels are reversed in 
the convolution in comparison to the beam labels. The parton distributions 
are passed from PYSTFU to PYSIGH through the array XPQ(KFL,KFLD) (see 
below) and stored in the array XSFX (N , KFL , KFLD) , where KFLD denotes the 
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helicity. ISIG(N, I) contains the information about the Nth-line in the event 
listing. The new entries ISIG(N,4) and ISIG(N,5) store the helicities of the 
partons at the beam and target side respectively. ISIG(N,6) is reserved but 
not used at the moment. 

SUBROUTINE PYSTFU(KF,X,Q2,XPQ) 

Purpose: to evaluate the helicity dependent parton distributions for given 
flavour (KF), x (X), and Q 2 (Q2) according to the selected parametrizations. 
New parameters: MSTP(177), MSTP(178), MSTP(179), MSTP(180), 
NSUB(ISUB) 

Internal Polarization switch: IP0L=MSTP(180) x (1-MSTP(179) ) 
Remarks: Call of PYSTFU returns x times the parton distribution functions 
for given flavour, x, and Q 2 for both helicities and an averaged (unpolar- 
ized) value. The values are stored in the array XPQ(KFL,KFLD) which has 
been enlarged from XPQ(-25:25) to XPQ(-25:25,0:2). Row XPQ(KFL,0) 
contains the unpolarized distributions, XPQ (KFL , 1) distributions correspond- 
ing to the positive helicity (relative to the hadron) and XPQ (KFL, 2) distri- 
butions for the negative helicity. The parton distributions are selected by 
switches described earlier (see Table p. The polarized distributions q±, g± 
are constructed from the unpolarized ones q, g (selected by old Pythia 
switches) and polarized parts Aq and Ag selected by MSTP(177) (see Table 
|l|). Four subroutines have been added to calculate the polarized distribu- 
tions. These are ALTSTI0, POLLO, and POLNLOQ The subroutines require data 
files altsti.dat, gestloa.dat, gestlob.dat, gestloc.dat, gestnloa.dat, 
gestnlob.dat, and gestnloc.dat. These files must be visible to FORTRAN 
open statement and therefore they have to be placed in an appropriate di- 
rectory. The files are supplied with the program. The parton distribution 
for fixed helicity are reconstructed from the polarized and unpolarized distri- 
butions according to q± — ^ (q ± Ag). When the polarized and unpolarized 
parts are combined together the program performs the unitarity check - if 
the resulting total distribution becomes negative for one helicity it is put to 
zero and the corresponding result for the other helicity is set to the value 

2 This subroutine has been written by G. Altarelli and J. Stirling. Used with permision 
from the authors. 

3 These subroutines have been written by T. Gehrmann and W.J. Stirling. Used with 
permision from the authors. 
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of the unpolarized part. Only polarized parametrizations for protons are 
implemented. Neutron parametrizations are obtained from them by isospin 
symmetry, the parametrizations for hyperons are constructed by naive SU(3). 
Charge conjugations is used to describe the corresponding antiparticles. 

The interface to the CERN parton distribution library was updated and 
works now with PDFLIB 4 and higher versions. The PAKPDF parton distribu- 
tion library is not supported anymore. Furhermore, several new unpolarized 
parton distribution sets have been included in Sphinx (see Table f|). 

When plotting ratios of polarized and unpolarized parton distribution func- 
tions one should have in mind that the distributions are often calculated from 
a two-dimensional grid by means of linear interpolation. These grids are 
usually logarithmic in x and Q 2 . Furthermore, the subroutines which give 
the polarized and unpolarized distributions may use different (x, Q 2 ) points 
and different interpolation routines. As a result, the asymmetry may not be 
smooth in the large-x region. Usually, this problem can be circumvented by 
using a parametrization of unpolarized parton distributions which is rather 
given by a fit with polynomials than by interpolation between data points. 
This method is usually slower but provides much smoother results. In gen- 
eral it is recommended to use the set of unpolarized parton distributions which 
corresponds to the polarized set. For the Gehrmann/ 'Stirling parametrizations 
these are the GRV94(LO) and MRS(AP)(NLO) sets, respectively. These sets 
are available in Sphinx as fits with polynomials by setting MSTP(51)=12 and 
MSTP(51)=13, respectively. 

SUBROUTINE PYSTPR(X,Q2,XPPR) 

Purpose: to calculate the unpolarized parton distributions of the proton; 
Remark: the implementation of the Cteq2 type unpolarized parton distri- 
butions was updated. The Interface to PYCTQ2 was taken from Pythia 5.7 
Several new unpolarized sets were included (see Table |). 

FUNCTION PYCTQ2 (Iset , Iprt ,X , Q) 

Purpose: to give the revised Cteq2 parton distribution sets with extended 
range in parametrized form. 

Remark: This function was taken from Pythia 5.7. 
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SUBROUTINE MRSAP (X , SCALE , UPV , DNV , USEA , DSEA , STR , CHM , BOT , GLU) 
Purpose: to calculate the unpolarized parton distribution functions set A' 
from Martin, Roberts and Stirling, x and Q are passed to the subroutine via 
the parameters X and SCALE. UPV, DNV, USEA, DSEA, STR, CHM, BOT, and GLU 
contain w va i, d val , u, d, s, c, b, and g repectively. 

Remarks: This subroutine was written by A.D. Martin, R.G. Roberts, and 
W.J. Stirling. For further information see ||. 



SUBROUTINE MRSR1 (X , Q2 , UPV , DNV , USEA , DSEA , STR , CHM , BOT , GLU) 
Purpose: to return the unpolarized parton distribution functions set Rl 
from Martin, Roberts and Stirling, x and Q 2 are passed to the subroutine 
via the parameters X and Q2. UPV, DNV, USEA, DSEA, STR, CHM, BOT, and GLU 
contain u Ya \, d va \, u, d, s, c, b, and g repectively. 

Remarks: This subroutine was written by A.D. Martin, R.G. Roberts, and 
W.J. Stirling. For further information see [JTOl . 



SUBROUTINE MRSR2 (X , Q2 , UPV , DNV , USEA , DSEA , STR , CHM , BOT , GLU) 
Purpose: to calculate the unpolarized parton distribution functions set R2 
from Martin, Roberts and Stirling, x and Q 2 are passed to the subroutine 
via the parameters X and Q2. UPV, DNV, USEA, DSEA, STR, CHM, BOT, and GLU 
contain u V3 \, d va \, u, d, s, c, b and g repectively. 

Remarks: This subroutine was written by A.D. Martin, R.G. Roberts, and 
W.J. Stirling. For further information see [FLO] . 



SUBROUTINE GRV94 (X , Q2 , UV , DV , DEL , UDB , SB , GL) 

Purpose: to return the leading order (1994) unpolarized parton distribution 
functions from Gliick, Reya, Vogt. x and Q 2 are passed to the subroutine 
via the parameters X and Q2. UV, DV, DEL, UDB, SB, and GL contain u Ya \, d va \, 
d — u, d + u, s, and g respectively. The heavy quark distributions c, b, and t 
are zero in the whole Q 2 -range. 

Remarks: This subroutine was written by M. Gliick, E. Reya, and A. Vogt. 



For further information see 111 



SUBROUT I NE DPL I ST ( ML I ST ) 

Purpose: to display the polarizations of the particle in the event listing. 
New parameters: KD(I) 
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Remarks: DPLIST is a modification of the Jetset subroutine LULIST. It is 
changed to display the polarization in the final listing. The sign displayed 
just behind the particle code denotes polarization with respect to the z-axis. 
When the sign is missing the particle has been treated as unpolarized. The 
information is taken from the vector KD(I) and transformed accordingly to 
('0V1')'2')— ','+','— '). The following format is chosen (polarization for the 
colliding hadrons, helicity for the partons resp.): 

• = ' ': no polarization/helicity 

• = '+': positive polarization/helicity 

• = ' — ': negative polarization/helicity 
SUBROUTINE DPEDIT(MEDIT) 

Purpose: to compress the vector KD(I), containing the polarization infor- 
mation, properly. 
New parameters: KD(I) 

Remarks: DPEDIT is a modification of the Jetset subroutine LUEDIT. 

SUBROUTINE ALTSTI (X , 02 , UPV , DNV , SEA , STR , CHM , BOT , TOP , GLU) 
Purpose: to return x times the polarized parton distributions evaluated at 
given x (X) and Q 2 (Q2) according to the parametrization of Altarelli&Stirling 
[||]. UPV denotes the valence distribution of up-quarks iA« val (i, Q 2 ), DNV for 
down-quarks xAd val (x, Q 2 ). SEA signifies the sea distribution xAq sca (x,Q 2 ). 
STR, CHM, BOT, and TOP label the distributions for the strange-, charm-, 
bottom-, and top-quark xAq(x, Q 2 ), q = s,c,b,t respectively. Finally GLU 
marks the gluon distribution xAg(x,Q 2 ). 

Remarks: ALTSTI requires the data file altsti .dat which has to be placed 
in an appropriate directory, altsti.dat is supplied with this program. 

SUBROUTINE P0LL0 ( IFLAG , X , Q2 , UVAL , DVAL , GLUE , QBAR , STR) 

Purpose: to return x times the polarized parton distributions evaluated 
at given x (X) and Q 2 (Q2) according to the parametrization of Gehrmann 
and Stirling lo set a (IFLAG=0), lo set b (IFLAG=1) or lo set c (IFLAG=2) g]. 
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UVAL, DVAL, GLUE, QBAR, STR contain Aw va i, Ad va i, A#, Ad = Am, and 5s 
repectively. 

Remarks: POLLO requires the subroutines RDARRY, POLINI and the data files 
gestloa.dat, gestlob.dat, gestloc.dat. Like this subroutine they were 
written by T. Gehrmann and W.J. Stirling and are used here with permission 
of the authors. For further informations about this package see 0. The data 
files have to be placed in the directory specified in POLINI and are supplied 
with this program. 

SUBROUTINE POLNLO (IFLAG , X , Q2 , UVAL , DVAL , GLUE , UBAR , DBAR , STR) 

Purpose: to return x times the polarized parton distributions evaluated at 
given x (X) and Q 2 (Q2) according to the parametrization of Gehrmann and 
Stirling nlo set a (IFLAG=0), nlo set b (IFLAG=l) or nlo set c (IFLAG=2) [|]. 
UVAL, DVAL, GLUE, UBAR, DBAR, STR contain Au vab Ad val , Ag, Au, Ad, and 
5s repectively. 

Remarks: POLNLO requires the subroutines RDARNLO, NLOINI and the data 
files gestnloa.dat, gestnlob.dat, gestnloc.dat. Like this subroutine 
they were written by T. Gehrmann and W.J. Stirling and are used here 
with permission of the authors. For further informations about this package 
see 0. The data files have to be placed in the directory specified in NLOINI 
and are supplied with this program. 

SUBROUTINE PDFSET (PARM , VALUE) 

Purpose: new dummy routine for new Pdflib interface; has to be removed 
if Pdflib is used. 

Remarks: This subroutine was taken from Pythia 5.7. 

SUBROUTINE STRUCTM (XX , QQ , UPV , DNV , USEA , DSEA , STR , CHM , BOT , TOP , GLU) 

Purpose: new dummy routine for new Pdflib interface; has to be removed 
if Pdflib is used. Remarks: This subroutine was taken from Pythia 5.7. 

4 Examples 

In the following we give two examples of a main program for a simulation 
with Sphinx and show the corresponding results. We considered longitudi- 
nal polarized proton-proton scattering in the CMS at y/s = 200 GeV and 
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selected the process qg — > qg for the partonic interaction. In the first ex- 
ample both beam and target are polarized in -^-direction, in the second 
example the target spin is reversed. With regard to the event listings is has 
to be mentioned that the displayed format differs from the real because we 
removed a few columns such that it fits in this text. In addition the event 
listing is cut after the hard interaction (denoted by • • •), i. e. behind line 8. 
In the omitted part there is no polarization information and it has the same 
format as the original Pythia listing. 

The information about the polarization flow is displayed as explained 
above right behind the flavour code KF. The first row contains the information 
about the beam particle. In addition to the Pythia labels the sign '+' 
behind the flavour code for the proton KF=2212 denotes the polarization in 
positive z-direction. Accordingly the sign '+' (' — ') in the second row signifies 
the polarization of the target in positive (negative) z-direction in the first 
(second) example. The third and fourth line represents the initiators of the 
initial state shower. In both examples they are an anti-d-quark and an u- 
quark with positive helicities for both partons. During the initial state shower 
the anti-quark becomes an gluon with negative helicity, whereas the target 
side parton remains a positive helicity u-quark. These partons undergo the 
hard interaction the resulting partons of which are displayed in the lines 
seven and eight. In the following the final state interaction takes place, 
i. e. the outgoing partons fragment, the unstable produced hadrons decay, 
etc. as long as only stable particles exist. In the final state polarization 
is not traced and consequently there is no polarization informations about 
these lines provided. This part of the listing is then again the same as the 
corresponding PYTHIA listing. 

4.1 The Main Programs 

4.1.1 First Example - parallel polarization 

C Example of a Main Program for event generating 

C in longitudinal polarized proton-proton-scattering 



C 
C 
C 
C 



E.g. : polarized p(+)p(+)-scattering in CMS 
at sqrt(s)=200 GeV 



24 



C This program has to be linked with 

C the programs SPHINX and JETSET7 . 3 , 

C the data files ALTSTI.DAT and R0SR0A.DAT, 

C and the CERN-Libraries . 

PROGRAM EXAMPLE 
C COMMON BLOCKS of SPHINX for event generation 



C0MM0N/LUDAT1/MSTU(200) ,PARU(200) ,MSTJ(200) , PAR J (200) 
COMMON/LUJETS/N , K (4000 , 5) , P (4000 , 5) , V (4000 , 5) 
C0MM0N/PYSUBS/MSEL,NSEL,MSUB(200) ,NSUB(200) , 
& KFIN(2, -40:40) ,CKIN(200) 

C0MM0N/PYPARS/MSTP(200) ,PARP(200) ,MSTI(200) , PARI (200) 
C0MM0N/PYINT5/NGEN (0 : 200 , 3) , XSEC (0 : 200 , 3) 

c============================================================: 

C Polarization Set-Up 

c============================================================: 

C Polarized simulation 

MSTP(180)=1 

C Beam positive polarized 

MSTP(171)=1 

C Target positive polarized 

MSTP(172)=1 

C Polarized parton distributions a la Gehrmann Stirling 

MSTP(177)=3 

C Polarized Initial State Shower 

MSTP(176)=1 
C with polarized parton distributions 

MSTP(175)=1 
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c========================================================== 

c========================================================== 

C Event Set-Up 

c========================================================== 

C Choice of parton processes "a la carte'' 

MSEL=0 

C Choice of polarized parton processes "a la carte'' 

NSEL=0 

C Choice of process: qg — > qg 

MSUB(28)=1 

C Process 28 polarized 

NSUB(28)=1 

C Number of generated events 

NEVENT=2000 

C kinematical cuts 

C P_T-cut (minimum) 

CKIN(3)=3. 
C___ P_T-cut (maximum) 

CKIN(4)=25. 

c========================================================== 

c========================================================== 

C Start of event generation 

c========================================================== 

C Initialisation 

CALL PYINITCCMS' , 'p' , 'p' ,200.) 
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Loop over events 
DO 100 I=1,NEVENT 



C Event generation 

CALL PYEVNT 

C List the first event with spininf ormation 

IF(I.EQ.l) CALL DPLIST(l) 

100 CONTINUE 

C Print cross section and histogram 

CALL PYSTAT(l) 

END 

4.1.2 Second Example — antiparallel polarization 

The second example is constructed by replacing MSTP(172)=1 by MSTP(172)=2 
in the first example, i. e. by switching the spin of the target. 

4.2 The Event Listings 

4.2.1 First example — parallel polarization 

SPHINX version 1.1 
** Last date of change: 18 Nov 1996 ** 

The Lund Monte Carlo - PYTHIA version 5.6 
** Last date of change: 3 Apr 1992 ** 

The Lund Monte Carlo - JETSET version 7.3 
** Last date of change: 10 Mar 1992 ** 

1******* PYINIT: initialization of PYTHIA routines ********** 
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I I 

I PYTHIA will be initialized for a p on p collider I 

I at 200.000 GeV center-of-mass energy I 

I I 



PYMAXI : summary of differential cross-section maximum search 



I ISUB Subprocess name I Maximum value I 



I II 

I28f+g->f+g I 4.7736E+00 I 

I 96 Semihard QCD 2 -> 2 I 1.9598E+02 I 

I II 



************** PYINIT: initialization completed ************* 

Event listing (summary) 



I 


particle/ jet KF 


P- 


.X 


P- 


-y 


p_z 


E 




m 




1 


!p+! 


2212+ 





000 





000 


99 


996 


100 


000 





938 


2 


!p+! 


2212+ 





000 





000 


-99 


996 


100 


000 





938 


3 


!d~! 


-1+ 


-o 


279 


-o 


107 


12 


583 


12 


587 





000 


4 


!u! 


2+ 


-0 


132 





469 


-3 


602 


3 


634 





000 


5 


!g! 


21- 


-2 


214 


-1 


299 


7 


340 


7 


776 





000 


6 


!u! 


2+ 


-0 


096 





342 


-2 


626 


2 


649 





000 


7 


!g! 


21 


-5 


168 


-1 


523 


5 


182 


7 


476 





000 


8 


!u! 


2 


2 


857 





566 


-0 


467 


2 


950 





006 
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1* PYSTAT: Statistics on Number of Events and Cross-sections * 



I III 

I Subprocess I Number of pointsl Sigma I 

I III 

I 1 1 (mb) I 

I III 

I N:o Type I Generated Tried I I 

I III 



I III 

I All included subprocesses I 2000 13360 I 7.385E-01 I 

I 28 f + g -> f + g I 2000 13360 I 7.385E-01 I 

I III 



* Fraction of events that fail fragmentation cuts = 0.00000 * 

4.2.2 Second example — antiparallel polarization 

SPHINX version 1.1 
** Last date of change: 18 Nov 1996 ** 

The Lund Monte Carlo - PYTHIA version 5.6 
** Last date of change: 3 Apr 1992 ** 

The Lund Monte Carlo - JETSET version 7.3 
** Last date of change: 10 Mar 1992 ** 

1********* PYINIT: initialization of PYTHIA routines ****** 
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I I 

I PYTHIA will be initialized for a p on p collider I 

I at 200.000 GeV center-of-mass energy I 

I I 



PYMAXI : summary of differential cross-section maximum search 



I II 
I ISUB Subprocess name I Maximum value I 

I II 



I II 

I28f + g->f + g I 4.7851E+00 I 

I 96 Semihard QCD 2 -> 2 I 1.9598E+02 I 

I II 



************ PYINIT: initialization completed ************** 



Event listing (summary) 



I 


particle/ jet KF 


P- 


_x 


P- 


-y 


p_z 


E 




m 




1 


!p+! 


2212+ 





000 





000 


99 


996 


100 


000 





938 


2 


!p+! 


2212- 





000 





000 


-99 


996 


100 


000 





938 


3 


!d~! 


-1+ 


-o 


279 


-o 


107 


12 


583 


12 


587 





000 


4 


!u! 


2+ 


-o 


132 





469 


-3 


602 


3 


634 





000 


5 


!g! 


21- 


-2 


214 


-1 


299 


7 


340 


7 


776 





000 


6 


!u! 


2+ 


-0 


096 





342 


-2 


626 


2 


649 





000 


7 


!g! 


21 


-5 


168 


-1 


523 


5 


182 


7 


476 





000 


8 


!u! 


2 


2 


857 





566 


-0 


467 


2 


950 





006 
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1* PYSTAT: Statistics on Number of Events and Cross-sections * 



I III 

I Subprocess I Number of pointsi Sigma I 

I III 

I 1 1 (mb) I 

I III 

I N:o Type I Generated Tried I I 

I III 



I III 

I All included subprocesses I 2000 12655 I 7.361E-01 I 

I 28 f + g -> f + g I 2000 12655 I 7.361E-01 I 

I III 



* Fraction of events that fail fragmentation cuts = 0.00000 * 
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